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ABSTRACT 

Context.  Generally,  there  are  tv\-o  procedures  lor  .solar  cycle  predictions:  the  empirical  methods  -  statistical  methods  based  on 
exirapolaiions  and  precursoi  methods  -  and  metlKxls  based  on  dynamo  models. 

Aims.  The  goal  of  the  present  analysis  is  to  forecast  the  strength  and  epochs  of  the  next  solar  cycle,  to  invesiigate  proxies  for  grand 
solar  minima  and  to  reconstruct  the  relaiive  sunsjxx  number  in  the  Maunder  mininuim. 

Methods.  We  calculate  the  asyntmetry  of  the  ascending  and  descending  solar  cycle  phases  (Method  1 )  and  use  this  parameter  as  a 
pi'oxy  for  solar  activity  on  longer  time  scales.  Further,  vs  e  correlate  the  relative  sunspot  numbers  in  the  ejXKhs  of  solar  activity  minima 
and  maxima  (Method  2)  and  estimate  the  parameters  of  an  autoregressive  moving  average  model  (ARMA,  Method  3).  Idnally,  the 
power  spectrum  of  data  obtained  with  the  Method  1  is  analysed  and  the  Meihixls  1  and  3  are  combined. 

Results.  Signatures  of  the  Maunder.  Dalton  and  Gleissberg  minima  were  found  with  Method  I .  A  period  of  about  70  years,  somewhat 
shorter  than  the  Glcissl>crg  period  was  identified  in  the  asymmetry  data.  The  maximal  smoothed  monthly  sunspot  number  during  the 
Maunder  minimum  was  reconstructed  and  found  to  be  in  the  range  0-35  (Method  1).  The  estimated  Wolf  number  (also  called  the 
relative  sunspot  number)  of  the  next  solar  maximum  is  in  the  range  88-102  (Metlxxl  2).  Method  3  predicts  the  next  solar  maximum 
between  201 1  and  2012  and  the  next  solar  minimum  for  2017.  Also,  it  forecasts  the  relative  sunspot  number  in  the  next  maximum  to 
he  90  ±  27.  A  combination  of  the  Methods  1  and  3  gives  for  the  next  solar  maximum  relative  sunspot  numhers  between  78  and  99. 
Conclusions.  The  asymmetry  parameter  provided  by  Method  1  is  a  good  proxy  for  solar  activity  in  the  past,  also  in  the  periods  for 
which  no  relative  sunspot  numbers  are  available.  Our  prediction  for  the  next  solai  cycle  No.  24  is  that  it  will  be  weaker  than  the  last 
cycle,  No.  23.  This  prediction  is  based  on  various  independent  mellxxls. 

Key  words.  Sun:  activity 


1.  Introduction 

Understanding  the  .solar  activity  cycle  rcmnins  a  key  unsolved 
problem  in  solar  physics  (along  with,  c.g.,  heating  of  the  solar 
corona  and  solar  Hares).  It  is  not  only  an  outstanding  theoretical 
problem,  hut  also  an  important  practical  issue,  since  solar  activ¬ 
ity  and  the  radiation  output  related  to  it  influences  the  biosphere, 
space  weather  and  technology  on  the  Earth  (Eddy  1 976:  Hoyt  & 
Schatten  1 997;  Feynman  &  Gabriel  2000;  Lang  2()()();  Soon  & 
Yaskell  2003;  Feynman  2007;  Hanslmeier  2007). 

At  present  it  is  still  not  fully  clear  whether  solar  activity 
is  purely  stochastic  or  weakly  chaotic  (Sti.x  2002;  Osscndrijver 
2003;  Rudiger  &  Hollcrbaeh  2004).  A  distinction  between  these 
two  charactens  (Carbonell  ct  al.  1 994;  Osscndrijver  &  Hoyng 
1 996;  Duhau  2003;  Charbonneau  2005;  Letcllicr  ct  al.  2006; 
Usoskin  ct  al.  2007;  Ycates  ct  al.  2008)  can  be  important 
for  the  interpretation  of  the  long-term  modulation  of  the  solar 
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activity  (e.g.,  the  Gleissberg  cycle  with  a  period  of  about 
lOO  years,  Feynman  &  Fougere  19S4;  Feynman  &.  Gabriel 
1 990;  Duhau  2003)  and  understanding  the  prolonged  periods 
of  extremely  low  activity,  as  in  the  Maunder  minimum  (e.g., 
Ruzmaikin  1 985;  Charbonneau  2005;  Hircmath  2006;  Volobucv 
2006).  In  this  context  a  very  important  but  difficult  task  is  to  pre¬ 
dict  a  possible  return  of  the  Maimder-type  grand  solar  minima. 

Solar  cycle  prediction  is  not  an  easy  task  and  there  exists  a 
variety  of  methods  based  on  different  observational  data  and  the¬ 
oretical  assumptions  (Wilson  1 994;  Beck  ct  al.  1 995;  Hathaway 
ct  al  1 999;  Osscndrijver  2003;  Schiissjer  2007).  Currently  there 
are  two  basic  classes  of  methods  for  the  .solar  cycle  predic- 
tion.s;  these  are  the  empirical  methods  and  methods  based  on 
dynamo  models  (c.g.,  Schiisslcr  2007).  The  empirical  methods 
(c.g.,  Hathaway  ct  al.  1 999)  can  be  further  devided  into  two 
subgroups,  the  statistical  methods  based  on  extrapolation  (e.g., 
Kane  20()7b)  and  precursor  methods  (e.g.,  Wilson  I99()b).  The 
other  class  of  methods  is  based  on  various  dynamo  models  (c.g., 
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Cameron  &  Seluisslcr  2007;  Choudhuri  ct  al.  2007;  Dikpali 
2007;  Jiang  et  al.  2007),  hut  can  al.so  be  combined  with  some 
precursor  features  (e.g.  polar  magnetic  field  of  the  Sun  around 
solar  minimum,  Svalgaard  ei  al.  2005).  Further,  various  com¬ 
binations  of  different  methods  were  proposed  (e.g.,  Hanslmeier 
ct  al.  1999;  Lanto.s*  2006a, b)  and  finally,  a  nonuni  form  latitude 
distribution  of  solar  activity  and  North-South  asymmetry  were 
also  used  for  predicting  future  solar  activity  (e.g.,  Javaraiah 
2007;  Kane  2007a). 

We  should,  however,  point  out  that  there  are  important  dif¬ 
ficulties  in  predieting  the  solar  cycle  using  mean-held  dynamo 
models  (Bushby  &  Tobias  2007;  Cameron  &  Schiissler  2008). 
The.sc  difficulties  arise  from  the  significant  modulation  of  the 
solar  activity  cycle  eau.sed  by  either  stochastic  or  deterministic 
processes.  Long-term  prediction  is  impossible  even  if  a  model 
used  for  the  prediction  would  be  correct  in  all  details.  Moreover, 
even  short-term  prediction  from  mean- field  models  is  problem¬ 
atic  because  of  fundamental  uncertainties  in  the  form  and  ampli¬ 
tude  of  the  transport  coefficient  and  nonlinear  respon.se.  All  this 
makes  the  fine  tuning  of  various  dynamo  models  controversial 
for  predicting  future  solar  activity. 

Generally,  it  can  be  said  that  for  the  next  solar  cycle  (No.  24) 
the  empirical  methods  predict  lower  amplitudes  (Sehatten  2005; 
Svalgaard  et  al.  2005;  Du  &  Du  2006;  Javaraiah  2007;  Hiremath 
2008)  and  the  methods  based  on  dynamo  models  higher  am¬ 
plitudes  (Dikpati  &  Gilman  2006;  Dikpati  et  al.  2006,  2008b). 
However,  this  distinction  is  not  a  strict  one,  and  in  both  classes 
of  methods  opposite  examples  can  be  found.  So  an  empirical 
method  based  on  geomagnetic  precursors  forecasts  a  high  am¬ 
plitude  for  the  next  solar  cycle  (Hathaway  &  Wilson  2006)  and  a 
solar  dynamo  model  predicts  a  low  amplitude  (Choudhuri  et  al. 
2007;  Jiang  et  al.  2007).  Finally,  various  combined  methods  pre¬ 
dict  a  low  amplitude  for  the  next  solar  maximum  (Sehatten  2003; 
Sehatten  Sc  Tobiska  2003;  Lantos  2()06a).  In  summary,  the  topic 
of  solar  cycle  prediction  is  still  inconclusive  to  a  large  extent, 
which  represents  additional  motivation  for  further  research. 

In  the  present  study  three  different  empirical  methods  for  so¬ 
lar  cycle  prediction  and  reconstruction  are  used.  We  calculate  the 
asymmetry  of  the  ascending  and  descending  solar  cycle  phases 
(Method  1)  and  use  this  parameter  as  a  proxy  for  solar  activity 
on  longer  time  scales.  Further,  we  correlate  the  relative  sunspot 
numbers  in  the  epochs  of  solar  activity  minima  and  maxima 
(Method  2).  Finally,  we  estimate  the  parameters  of  an  autore¬ 
gressive  moving  average  model  (ARMA,  Method  3). 

With  Method  1  we  search  for  indications  of  grand  solar  min¬ 
ima  in  the  past.  It  also  enables  the  reconstruction  of  the  relative 
sunspot  number  in  the  Maunder  minimum.  With  ihis  method  a 
prediction  for  the  amplitude  of  the  next  solar  cycle  is  possible 
only  if  the  epochs  of  future  minima  and  maxima  can  be  inde¬ 
pendently  estimated.  Method  2  provides  a  prediction  of  the  am¬ 
plitude  of  the  next  activity  maximum  on  the  basis  of  the  rela¬ 
tive  sunspot  number  in  the  preceding  minimum.  Finally,  with 
Method  3  predictions  of  the  amplitude  and  epochs  of  the  next 
solar  cycle  are  possible. 

2.  The  data  set 

Our  ba.sie  data  set  consists  of  the  epochs  of  .solar  minima  and 
maxima,  7 min  and  /’max.  respectively,  from  the  year  1610  up  to 
now.  The  starting  point  was  the  table  published  by  Gleissberg 
ct  al.  (1979).  The  epochs  of  solar  minima  and  maxima  were 
determined  by  Wolf  for  the  years  up  to  1889,  by  Wolfer  for 
the  years  up  to  1923,  by  Brunner  for  the  years  up  to  1944, 
and  by  Waldmcicr  for  later  years.  We  have  cheeked,  slightly 


corrected  and  enlarged  this  data  set  using  monthly  smoothed 
relative  sunspot  number  beginning  from  1750.  This  data  .scries 
ean  be  found  at  the  solar  influences  data  analysis  center  (SfDC) 
of  the  Royal  Observatory  of  Belgium  (SlDC-tcam).  The  data  set 
used  in  the  present  work  is  reproduced  in  Table  1  (we  note  that 
the  asymmetries.  A,  A  t ,  A2,  and  A3,  written  in  the  right-hand  part 
of  the  table,  are  calculated  according  to  Method  1,  described  in 
the  next  section). 

Wc  now  hriefly  discuss  the  errors  in  the  determination  of  the 
epochs  of  solar  minima  and  maxima,  7  min  ^nd  7',nax*  According 
to  Gleissberg  ct  al.  (1979)  these  errors  are  in  modern  times  of 
the  order  of  0.2  years  and  for  the  Maunder  minimum  about  three 
times  larger,  i.c.  slightly  more  than  half  a  year.  In  agreement  with 
that.  Wolf  himself  (1893)  estimated  the  errors  of  the  minimum 
and  maximum  epochs  in  the  17th  century  to  be  less  than  1  year. 

3.  The  reduction  methods 

3.1.  Method  1:  the  asymmetry  A  of  the  ascending  (T) 
and  descending  (U)  solar  cycle  phases 

The  a.symmetry  of  the  ascending  and  the  descending  .solar  cycle 
phases  is  determined  according  to  Gleissberg  ct  al.  ( 1 979).  Using 
the  time  of  a.seent  from  the  minimum  epoch  to  the  maximum 
epoch  (T)  and  the  time  of  descent  from  the  maximum  epoch  to 
the  next  minimum  epoch  (U)  the  asymmetry  is  calculated  with; 


The  asymmetry  A  is  positive  if  the  activity  maximum  occurs  in 
the  first  half  of  the  cycle  and  negative  if  it  falls  in  the  second 
half.  The  A  values  were  calculated  using  Eq.  ( I )  for  all  available 
data  and  listed  in  Table  I . 

Further,  various  smoothing  procedures  were  applied; 


/!(/- I)  + 2/1(0 +  /^(/+ 1) 

/l,(0=  ^ 

(2) 

^  A(i  -  \)  +  AO)  +  AO  +  \) 

MO)  =  ^ 

(3) 

a.s  well  as  the  smoothing  method  of  Gleissberg  et  al.  (1979). 
With  this  method  the  moving  averages  of  every  four  consecu¬ 
tive  A  values  are  computed  and  from  this  series  again  the  means 
of  every  two  adjacent  values  arc  calculated.  In  such  a  way  the 
smoothed  values  A3  were  obtained: 

,  AiiiO+AnU) 

MO)  =  - ^ - .  (4) 

with 

^  A(i-2)+Aii-])+A0)  +  A0+^) 

A31  (/)  = - ^ ,  (5) 

4 

and 

X  /•^  A(i  -  1)  +  A(/) -I- A(i -I-  1)  +  A(/ -I- 2) 

A32(0  =  - - - (6) 

4 

In  Eqs.  (2)-(6),  A(/  -  1 )  is  the  asymmetry  of  the  solar  cycle  pre¬ 
ceding  to  the  cycle  i  and  A(/  -h  1 )  the  asymmetry  of  the  solar  cycle 
following  the  nh  cycle.  In  a  similar  way  A(/  2)  and  A(/  +  2)  are 
also  defined. 

All  these  asymmetry  values.  A,  Ai,  A2»  and  A3,  calculated  in 
the  above  de.scribed  way  are  presented  in  the  right-hand  part  of 
Table  1.  From  the  smoothing  procedures  it  is  clear  that  the  first 
and  the  last  entries  cannot  be  calculated  for  Ai  and  A2,  which 
holds  for  the  first  two  and  last  two  entries  in  the  case  of  A3. 
“NaN"  stands  for  no  data  available  or  the  value  cannot  be  cal¬ 
culated  and  the  epoch  of  the  actual  solar  minimum  is  only  an 
approximation. 
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Table  I.  The  solar  cycle  number,  the  eptxhs  of  solar  minima  and  maxima  (7‘min  and  respeciively),  the  extreme  values  of  the  monthly  smtK)thed 
lelatue  sunspot  numbers  in  corresponding  epcKhs  and  respectively)  and  ihe  asymmetry  values  A|,  A2,  and  A.i). 


Cycle  No. 

Tnitn  (year) 

T'max  (year) 

^min 

/?n.ax 

A 

A. 

A2 

Al 

-12 

1610.8 

1615.5 

NaN 

NaN 

-0.1460 

NaN 

NaN 

NaN 

-11 

1619.0 

1 626.0 

NaN 

NaN 

0.0670 

-0.0030 

-0.0263 

NaN 

-10 

1634.0 

1639.5 

NaN 

NaN 

0.0000 

0.0668 

0.0890 

0.0599 

-9 

164.‘).0 

lf)49.0 

NaN 

NaN 

0.2000 

0.1227 

0.()970 

0.0395 

-8 

1655.0 

1 660.0 

NaN 

NaN 

0.0910 

0.0122 

-0,0140 

-0.0230 

-7 

1666.0 

1675.0 

NaN 

NaN 

-0.3330 

-0.1688 

-0.1140 

-0.0385 

-6 

1679.5 

1685.0 

NaN 

NaN 

-0.1000 

-0.0893 

-0.0857 

-0.0618 

-5 

1689.5 

1693  0 

NaN 

NaN 

0.1760 

0.0452 

0.0017 

-0.0501 

-4 

1698.0 

1705.5 

NaN 

NaN 

-0.0710 

-0.0110 

0.0090 

0.0240 

-3 

1712.0 

1718.2 

NaN 

NaN 

-0.0780 

0.0028 

0.0297 

0.0624 

-2 

1723.5 

1727.5 

NaN 

NaN 

0.2380 

0. 1 357 

0.1017 

0.0625 

-I 

1734.0 

1738.7 

NaN 

NaN 

0.1450 

0. 1 223 

0.1147 

0.0606 

0 

1745.0 

1750.3 

NaN 

92.6 

-0.0390 

-0.0145 

-0.0063 

0.0594 

1 

1755.2 

1761.5 

8.4 

86.5 

-0.1250 

-0.0035 

0.0370 

0.0921 

2 

1766.4 

1769.7 

112 

115.8 

0.2750 

0.1988 

0.1733 

0.1876 

3 

1775.5 

1778.4 

7.2 

158.5 

0.3700 

0.3788 

0.3817 

0.2562 

4 

1784.7 

1788.1 

9.5 

141.2 

0.5000 

0.3138 

0.2517 

0.2329 

5 

1798.3 

1805.1 

3.2 

49.2 

-0.1  150 

0.0870 

0.1543 

0.1314 

6 

1810.5 

1816.4 

0.0 

48.7 

0.0780 

-0.0510 

-0.0940 

0,031  1 

7 

1823.3 

1829.9 

0.1 

71.5 

-0.2450 

-0.0247 

0.0487 

0.0551 

8 

1833.9 

1837.2 

7.3 

146.9 

0.3130 

0.1613 

0.1  107 

0. 1 263 

9 

1843.5 

1848,1 

10.6 

131.9 

0.2640 

0.2772 

0.2817 

0.2336 

10 

1856.0 

1860.1 

3.2 

98.0 

0.2680 

0.3060 

0.3187 

0.2915 

11 

1867.2 

1870.6 

5.2 

140.3 

0.4240 

0.3058 

0.2663 

0.2779 

12 

1879.0 

1884.0 

2.2 

74.6 

0.1070 

0.2497 

0.2973 

0.2937 

13 

1890.2 

1894.0 

5.0 

87.9 

0..3610 

0.2818 

0.2553 

0  2680 

14 

1902.1 

1906.1 

2.7 

64.2 

0.2980 

0.2863 

0,2823 

0.2337 

15 

1913.5 

1917.6 

1.5 

105.4 

0.1880 

0.1857 

0.1850 

0.2224 

16 

1923.6 

1928  3 

5.6 

78.1 

0.0690 

0.1585 

0.1883 

0.2226 

17 

1933  7 

1937.3 

3.5 

119.2 

0.3080 

0.2595 

0.2433 

0.2381 

18 

1944.1 

1947.4 

7.7 

151.8 

0.3530 

0.3177 

0.3060 

0.2740 

19 

1954.3 

1958.2 

3.4 

201.3 

0.2570 

0.2885 

0.2990 

0.2988 

20 

1964.8 

1968.9 

9.6 

110  6 

0.2870 

0.2797 

0.2773 

0.3050 

21 

1976.3 

1980.0 

12.2 

164.5 

0.2880 

0.3215 

0.3327 

0.3187 

22 

1986  7 

1989.5 

12.3 

158.5 

0.4230 

0.3577 

0.3360 

NaN 

23 

1996.4 

2000.3 

8.0 

120.8 

0.2970 

NaN 

NaN 

NaN 

24 

2008.5 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

NaN 

3.2.  Method  2:  the  correlation  of  the  relative  sunspot 
numbers  in  the  solar  minima  Rfj,in  and  maxima 

I'his  method  is  based  on  the  assumed  linear  relationship  between 
relative  sunspot  numbers  in  the  minimum  (at  the  beginning  of 
the  activity  cycle,  /?miti)  ^ind  in  maximum  epochs  (/?,nax)  of  solar 
cycles  (c.g.,  Wilson  1990b).  For  the  correlation,  series  of  the  ex¬ 
treme  values  (measured  in  minimum  and  maximum  epochs)  of 
the  monthly  smoothed  relative  sunspot  number,  beginning  from 
1750  (Table  1),  were  used.  This  method  enables  predicting  of 
/?,„a.x  on  the  basis  of  the  preceding 

3.3.  Method  3:  the  autoregressive  moving  average  model, 
ARM  A 

F.specially  interesting  for  forecasting  are  models  where  the  true 
state  r(/)  and  its  forecasted  value  X{t)  deviate  from  each  other  by 
e  which  is  a  Gaussian  distributed  random  variable,  That  means 
that  the  dilference  .v(/)- Jc(0  between  true  state  and  the  predicted 
value  does  not  contain  information  on  the  stochastic  process  that 
produced  the  time  scries  jc(0.  Such  a  general  process  can  be  de- 
.scribed  by  the  equation 

x{i)  =  Yj  ~  ~ 

,.i  /-I 


with  e  IV.V(0,(7-^)  and  is  called  an  ARMA|/:>, (autoregres¬ 
sive  moving  average)  process  (Brockwell  &  Davies  1996).  An 
ARMA  process  is  a  superposition  of  an  autoregressive  process 
(AR[p])  and  a  moving  average  process  (MA[^]).  The  parame¬ 
ters  p  and  q  are  the  orders  of  the  process,  i.e.,  they  define  which 
time  step  values  from  the  past  influence  the  present  state.  The 
coefficients  and  hj  define  how  strong  this  inlluence  from  the 
past  is.  The  random  variables  e  obey  a  white  noise  (Gaussian) 
probability  distribution  WN{0,(t^)  with  the  maximum  at  0  and 
a  variance  of  rr^.  We  use  such  an  ARMA  process  to  model  the 
sunspot  number. 

In  the  following  we  describe  briefly  the  procedure  for  es¬ 
timating  the  ARMA  parameters  a,  and  bj  and  determining  the 
orders  p  and  q  from  the  lime  scries  of  sunspot  numbers.  More 
details  can  be  found  in  Durbin  &  Koopman  (2001)  and  Gardner 
et  al.  (1980).  First  we  use  initial  estimates  from  the  Yule- Walker 
equations,  i.c.,  we  assume  q  =  0.  Second,  the  residues  from  such 
an  ARIp]  fit  are  used  as  first  gues.scs  for  6(/).  Third,  we  minimi/c 

^(e(0)‘  =  ^  ( v(/)  -  [  Y  zLi  ~  -^^1} 

III  j 

to  obtain  belter  estimates  for  the  parameters  a,  and  h,.  We  de¬ 
cide  on  the  model  orders  p  and  q  by  increasing  the  orders  un¬ 
til  the  residues  arc  compatible  with  white  noise,  i.c.,  additional 
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Fig.  1.  The  asymmeiry  parameter  A  as  a  function  of  the  cycle  number. 
The  parameters  of  the  linear  lea.si-square  (it  are  also  given. 

parameters  do  not  add  additional  information  on  the  process  to 
the  model. 

Once  we  obtain  the  ARMA  model  on  the  basis  of  the  last 
358  yearly  values  we  use  it  for  forecasting  the  next  12  yearly 
values.  The  yearly  values  of  the  relative  sunspot  number  for 
the  period  1650-1699  were  taken  from  Eddy  (1976)  and  since 
the  year  1700  from  Waldmeier  (1961).  We  note  that  the  val¬ 
ues  published  by  Waldmeier  (1961)  were  reproduced  by  Hoyt  & 
Schalten  (1998,  the  Rz  values  in  their  Appendix  2).  The  eonlein- 
porary  data  are  taken  from  the  Solar  InHuenecs  Data  Analysis 
Center  (SIDC)  of  the  Royal  Ob.servatory  of  Belgium  (http:// 
sidc.oma.be). 

4.  Results 


Fig.  2.  Similar  lo  Fig.  1,  here  prescnied  for  the  smoothed  values  A^. 


4. 1.  Method  1:  the  asymmetry  A 

First  we  analyse  the  asymmetry  parameters  as  a  function  of  the 
cycle  number,  i.e.,  the  time.  These  dependencies  are  presented 
in  Figs.  I  and  2  for  the  asymmetry  parameters  A  and  A3,  re¬ 
spectively.  In  Fig.  3  the  power  .spectrum  as  a  function  of  the  in¬ 
verse  cycle  number,  i.e.  time,  for  the  asymmetry  data  series  from 
Fig.  1  is  presented.  The  Fourier  analysis  following  the  method  of 
Deeming  (1975)  was  applied  to  the  asymmeti7  parameter  A  as  a 
function  of  the  inverse  cycle  number.  The  highest  peak  in  Fig.  3 
corresponds  to  the  period  of  about  70  years,  which  is  slightly 
shorter  than  the  period  of  the  Glei.ssberg  cycle. 

The  asymmetry  parameter  A3,  smoothed  according  to  the 
method  of  Glei.ssberg  ct  al.  (1979),  is  presented  as  a  function 
of  the  solar  cycle  number  in  Fig.  2.  The  minima  in  the  smoothed 
curve  clearly  indicate  the  Maunder  minimum  (cycles  -9  to  -A; 
1645-1715),  the  Dalton  minimum  (cycles  5  lo  7;  1800-1830) 
and  the  modem  (Gleissbcrg)  minimum  (cycles  12  to  14;  1880- 
1915).  So,  the  a.symmctry  parameter  A3  is  a  good  proxy  for  grand 
solar  minima. 

Concerning  the  asymmetry  data  presented  in  Fig.  2  we 
should  make  two  eomments.  (i)  The  trend  of  increasing  asym¬ 
metry  with  lime  is  consistent  with  the  well-known  fact  that 
during  the  analy.scd  time  interval  the  solar  activity  also  in- 
crca.scd  on  average  (c.g.,  Schatten  2003,  Fig.  1 ;  McCracken  et  al. 
2004,  Fig.  le).  The  Gleissberg  and  the  Sehwabe  cycles  are  su¬ 
perimposed  on  this  increasing  trend  in  solar  activity.  We  ean 


Fig. 3.  The  power  spectrum  as  a  function  of  the  inverse  cycle  number 
for  the  data  scries  presented  in  Fig.  1 .  The  highest  peak  corresponds  to 
the  period  of  about  70  years. 

understand  the  increasing  trend  in  the  asymmetry  in  terms  of 
a  linear  relation  between  the  cycle  amplitude  and  correspond¬ 
ing  a.symmetry,  as  will  be  discussed  later  in  this  section,  (ii)  The 
minimum  of  the  asymmetry  curve  related  to  the  Gleissberg  min¬ 
imum  occured  in  cycles  14  to  16,  rather  than  in  cycles  12  to  1 4, 
when  the  Gleissberg  minimum  took  place.  For  this  time  shift  we 
have  no  simple  explanation,  except  that  there  are  some  indica¬ 
tions  that  the  Gleissberg  minimum  was  prolonged  for  one  or  two 
decades  (i.e.,  in  the  1920s  and  1930s). 

We  now  consider  the  four  asymmetry  parameters  as  a  func¬ 
tion  of  the  highest  relative  sunspot  number,  /?niax.  represented 
by  the  smoothed  monthly  values  (Table  1).  The  correlation 
coefficients  are  0.65,  0.64,  0.57,  0.53.  for  the  values  of  A,  Aj, 
A2,  and  A3,  respectively.  We  now  apply  the  linear  least-square  fit 
to  these  data  scries  of  the  general  form: 

A=a  +  /?/?max,  (9) 

where  a  and  h  are  the  fit  parameters.  Two  examples  of  these  de¬ 
pendencies  arc  presented  in  Figs.  4  and  5,  for  the  eases  At  (/?  max) 
and  A3(/?max)>  respectively.  We  now  .solve  Eq.  (9)  for 

-  M  -  a) 

^max 


h 


(10) 
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which  enables  a  reconstruction  of  the  relative  sunspot  number 
in  the  epochs  when  the  asymmetry  is  known.  For  the  Maunder 
minimum  we  gel  values  =  0,  /?nia.x  =  23,  and  /^,na\  =  35, 
for  diftcrenl  combinations  of  the  fit  parameters  a  and  h  and  the 
asymmetry  values  from  Table  1 ,  We  stress  that  in  some  cases  the 
reconstructed  relative  sunspot  numbers  are  negative,  so  that  not 
all  combinations  of  the  lit  parameters  and  asymmetry  values  arc 
possible.  It  is  also  important  to  note  that  these  arc  only  crude 
reconstructed  values  of  the  relative  sunspot  number.  The  recon- 
.siriicicd  values  for  the  Maunder  minimum  arc  lower  than  any 
peak  /?,nax  values  in  the  solar  cycles  thereafter,  as  can  be  seen  in 
Table  I . 


4.2.  Method  2:  the  correlation  of  Rmin  sind 

In  Fig.  6  we  present  a  relationship  between  /?,nin  and  /?,nax  for 
all  solar  cycles  w  ith  available  data  ('lablc  1 ).  As  can  be  seen  in 
the  ligure  the  solar  cycle  No.  19  represenis  a  kind  of  an  outlier, 
outside  of  the  bulk  of  data  points.  The  solar  activity  maximum 
in  1958  had  the  highest  /?n,ax  value  which  was  preceded  by  the 
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Fig,  6,  The  peak  smoothed  monthly  rchilivc  sunspot  number  as  a  func 
tion  of  the  same  quantity  in  the  preceding  .solar  minimum,  for  all  solar 
cycles  with  the  data  available.  The  parameters  ol  the  linear  least-square 
fit  are  also  given. 


unusually  low  minimum.  So,  in  the  further  analysis  we  exclude 
the  data  point  from  the  solar  cycle  No.  1 9.  That  solar  cycle  w  as  a 
special  case,  possibly  even  an  anomaly  (Wilson  1990a),  not  only 
because  of  the  extremely  low  minimum  at  its  beginning,  but  also 
due  to  extreme  North-South  asymmetry  of  the  activity,  as  re¬ 
ported  by  Temmcr  cl  al.  (2006).  This  pronounced  North-South 
asymmetry  might  be  a  signature  of  a  weak  mixed-parity  dynamo 
mode  (Sokoloff  &  Nesme-Ribes  1994).  One  mixed-parity  solu¬ 
tion  is  related  to  ordinary  Schwabe  cycles,  while  the  other  one 
is  relevant  for  grand  minima.  It  is  interesting,  however,  that  the 
anomaloLi.s  behaviour  of  cycle  No,  19  cannot  be  recognized  in 
the  cycle  asymmetry  data  (Figs.  1  and  2). 

Without  the  data  point  from  .solar  cycle  No.  19  w'e  obtain  the 
following  linear  relationship: 

/^nm  =67.4(±l().6)+6.9(±1.5)/?ni,n,  (H) 

w  ith  the  correlation  coefficient  of  0.72.  Equation  (11)  enables  a 
prediction  of  the  amplitude  of  the  next  solar  activity  maximum, 
if  the  /?niin  i''  known  (or  can  be  estimated).  For  the  estimated  /?,n,n 
values  in  the  current  solar  minimum  of  5,  4,  and  3,  w'c  obtain 
/?n,ax  values  of  1 02, 95,  and  88  for  the  next  maximum,  according 
to  Ecj.  (II). 

For  comparison,  w'e  repeal  one  of  the  formulae  obtained  by 
Wilson  ( 1 990b): 

/^max  =  81.7 -h  5.46 (12) 
4.3.  Method  3:  the  ARMA  model 

After  successful  fitting  of  a  model  {p  =  6,  ^  =  6)  the  estimated 
process  is  used  to  predict  the  key  quantities  for  the  next  solar 
cycle.  We  estimate  on  the  basis  of  the  yearly  predictions  that  the 
next  solar  activity  maximum  w'ill  take  place  between  201 1  and 
2012,  while  the  minimum  thereafter  should  occur  in  2017.  The 
model  predicts  the  height  of  the  next  maximum  in  the  sunspot 
number  to  be  /?,iiax  =  90  ±  27.  The  observed  yearly  values  of 
the  relati\e  sunspot  number  and  the  predicted  ones,  using  the 
ARMA  method,  arc  presented  in  Fig.  7. 
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Fig.  7.  The  observed  (marked  by  x  symbol  and  the  full  line)  and  pre¬ 
dicted  (marked  by  diamond  symbol  and  the  dotted  line)  yearly  values 
of  the  relative  sunspot  number  using  the  ARMA  method. 


4.4.  Combined  method  1  (the  asymmetry  A)  and  method  3 
(the  ARMA  model) 

F^inally,  we  can  combine  the  ARMA  ineihoci  (Method  3)  and  the 
asymmetry  method  (Method  I )  to  predict  the  strength  of  the  next 
solar  maximum  using  the  hypothesi.s  that  the  actual  solar  mini¬ 
mum  occurr.s  at  7’min  =  2008.5.  Prom  the  ARMA  model  we  get 
/  in:ix  =  20l2.0and  Z’min  =  20 17.0  (Sect.  4.3.).  Now  using  Eq.  (1) 
From  Sect.  3.1.,  the  a.symmetry  parameter  can  be  calculated  as 
A  =  0.176.  Inserting  this  value  into  Eq.  (10)  we  obtain  for  Z^niax 
the  values  99,  93,  92,  and  78,  for  the  A,  A|,  A?,  and  A?  cases, 
respectively. 

5.  Summary,  discussion  and  conclusions 

Signatures  of  the  Maunder  (1645-1715),  Dalton  (1800-1830) 
and  Gleissherg  (1880-1915)  minima  were  found  with  Method  1 . 
A  period  of  about  70  years,  somewhat  shorter  than  the 
Gleissherg  period  was  identified  in  the  power  spectrum  of  the 
asymmetry  data  as  a  function  of  time.  The  maximal  smoothed 
monthly  sunspot  number  during  the  Maunder  minimum  was  re¬ 
constructed  and  found  to  be  in  the  range  0-35  (Method  1 ).  The 
forcasted  strength  of  the  next  solar  maximum  is  in  the  range 
88-102  (Method  2).  The  ARMA  method  (Method  3)  predicts 
the  next  solar  maximum  between  201 1  and  2012  and  the  next 
solar  minimum  for  2017.  Also,  it  forecasts  the  relative  sunspot 
number  in  the  next  maximum  to  be  90  ±  27.  So,  we  forecast  a 
relatively  weak,  but  also  a  relatively  short  24ih  solar  cycle.  We 
stress,  however,  that  the  temporal  uncertainty  in  forecasting  7‘min 
and  with  ARMA  is  limited  hy  the  fact  that  we  used  here  the 
yearly  sunspot  number  values.  In  a  subsequent  paper  we  plan 
to  use  the  monthly  values  which  will  improve  the  precision  of 
the  predicted  epochs.  An  ARMA  proce.ss  is  a  linear  stationary 
stochastic  process.  The  sunspot  number  is  probably  the  result 
of  a  more  complex  system.  By  modelling  this  system  with  an 
ARMA  model  we  assume  that  this  system  is  only  weakly  non¬ 
stationary  and  non-linear.  A  longer  lime  .scries  of  the  sunspot 
numbers  is  requested  to  show  whether  this  assumption  is  cor¬ 
rect.  A  combination  of  Methods  I  and  3  gives  for  the  next  solar 
maximum  relative  sunspot  numbers  between  78  and  99. 

Our  prediction  for  the  next  solar  maximum  is  In  the  range 
from  /?max  =  78  to  /^max  =  102,  depending  on  the  chosen 
method,  combination  of  the  methods  and  precursor  parameters. 


Many  other  predictions  obtained  applying  the  empirical  or  com- 
billed  methods  are  within  this  span  (Schatten  2003,  2005:  Kane 
20()7b),  just  below  our  lower  limit  (Svalgaard  et  al.  2005; 
Javaraiah  2007),  slightly  above  our  higher  limit  (Lantos  2006a; 
Du  Si  Du  2006;  Hiremalh  2008),  or  significantly  above  it  (Kane 
2007a).  Our  prediction  is  in  qualitative  agreement  with  a  forecast 
based  on  a  dynamo  model  (Choudhuri  et  al.  2007;  Jiang  et  al. 
2007).  It  is  in  a  qualitative  and  quantitative  disagreement  with 
predictions  based  on  other  dynamo  models  (Dikpaii  &  Gilman 
2006;  Dikpati  et  al.  2006,  2008b)  and  with  one  empirical  fore¬ 
cast  (Hathaway  &  Wilson  2006). 

Concerning  the  start  of  the  cycle  No.  24,  we  can  now  ignore 
or  even  reject  some  earlier  prediciions  w'hich  put  It  in  early  2007 
(Du  2006;  Lantos  2006a).  Our  prediction  that  the  next  maximum 
should  occur  between  201 1  and  2012  is  roughly  in  agreement 
with  another  forecast  (Lantos  20()6a)  obtained  with  a  combined 
method. 

The  a.symmetry  parameter  provided  by  Method  1  is  a  good 
proxy  for  solar  activity  in  the  past,  also  in  the  periods  for  which 
no  relative  sunspot  numbers  are  available.  Our  prediction  for  the 
next  solar  cycle  No.  24  is  that  it  will  be  weaker  than  the  last 
cycle,  No.  23.  This  prediction  is  based  on  various  Independent 
methods  and  is  consistent  with  only  a  small  probability  that  the 
high  level  of  solar  activity  during  the  last  65  years  will  continue 
in  the  next  decades  (Schatten  &  Tobiska  2003;  Solanki  et  al. 
2004;  Usoskin  et  al.  2007).  On  the  other  hand,  we  should  men¬ 
tion  the  possibility  that  the  commonly  used  sunspot  number  se- 
rie.s  for  the  last  165  years  is  not  calibrated  correctly  (Giver  & 
Svalgaard  2007;  Svalgaard  &  Cliver  2007).  In  the  case  that  this 
will  prove  true,  the  secular  increase  of  the  solar  activity  in  the 
20th  century  would  be  called  into  question  and  our  prediction 
for  the  next  solar  maximum  would  be  affected  too.  However, 
the  times  of  activity  minimn  and  maxima  will  not  change.  We 
plan  to  include  the  analysis  of  the  corrected  sunspot  number  se¬ 
ries  in  our  future  research.  AI.so,  the  existing  long-term  series  of 
and  data  (e.g..  Beer  ct  al.  2006;  Vonmoos  el  al.  2006; 
Abreu  et  al.  2008)  might  be  suitable  for  an  application  of  the  em¬ 
pirical  methods.  It  is  not  surprising  that  in  the  ^^Be  time-series 
the  Maunder,  Dalton  and  Gleissherg  grand  minima  were  clearly 
identified  (McCracken  ct  al.  2004)  and  it  would  be  interesting 
to  sec  if  these  data  would  forecast  upcoming  grand  minima.  An 
analysis  of  correctly  calibrated  data  sets  from  these  terrestrial 
proxies  could  also  provide  a  better  insight  into  whether  or  not 
the  system  is  weakly  non-stationary  and  weakly  non-linear,  as 
earlier  di.scusscd. 

Further  we  can  raise  the  question  of  whether  the  cycle  asym¬ 
metry,  as  defined  by  Gleissherg  et  al.  (1979)  and  applied  in  the 
present  analysis  (Method  1),  is  at  all  an  important  proxy  of  the 
solar  activity  on  the  larger  temporal  scale.  We  propose  that  it  is 
indeed  the  case,  in  spite  of  some  other  opinions.  For  instance, 
Hathaway  ct  al.  ( 1994)  claimed  that  such  an  asymmetry  is  not 
only  not  very  important,  but  also  that  its  variation  from  cycle  to 
cycle  is  small  and  not  indicative  of  .solar  maxima.  In  the  present 
analysis  we  present  evidence  for  an  opposite  conclusion. 

A  final  question  is  related  to  the  possible  phy.sical  interpre¬ 
tation  (or  Justification)  of  the  study  of  the  asymmetry  param¬ 
eter.  Our  analysis  of  the  a.symmetry  and  its  proportionality  to 
the  cycle  strength  is  obviously  related  to,  but  is  not  the  same 
as  the  Waldmeier  effect.  This  etfect  states  that  stronger  cycles 
have  shorter  rising  times  (e.g.,  Wilson  199()c;  Hathaway  el  al. 
2002;  Lantos  2006b)  and  is  one  of  Waldmcicr’s  rules  (for  a  de¬ 
tailed  review  .sec  e.g.  Beck  ct  al.  1995).  The  Waldmeier  effect 
was  found  to  be  a  result  of  (or  at  least  consistent  with)  var¬ 
ious  dynamo  models,  starting  with  non-linear  dynamo  waves 
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(Sli\  1972).  various  flticlualioii.s  in  ihc  mean  field  dynamo  theory 
(Hoyng  1993;  Osscndrijvcr  &  Hoyng  1996)  up  lo  the  more  re¬ 
cent  fine  tuning  of  the  flux  transport  dynamo  models  (Cameron 
<&  Schiisslcr  2007).  On  the  other  hand,  it  was  even  considered 
that  the  Waldmcier  ctfecl  might  be  an  arlifacl  of  the  definition 
of  the  rehilive  sunspot  number  (Dikpali  cl  al,  20()8a),  since  it  is 
much  less  pronounced  in  the  sunpol  area  data  sets. 
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